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Preface 


Standard  Reference  Materials  (SRM's)  as  defined  by  the  National  Institute  of  Standards 
and  Technology  (NIST),  formerly  the  National  Bureau  of  Standards  (NBS),  are  well- 
characterized  materials,  produced  in  quantity  and  certified  for  one  or  more  physical  or 
chemical  properties.  They  are  used  to  assure  the  accuracy  and  compatibility  of  measure- 
ments throughout  the  Nation.  SRM's  are  widely  used  as  primary  standards  in  many  diverse 
fields  in  science,  industry,  and  technology,  both  within  the  United  States  and  throughout 
the  world.  They  are  also  used  extensively  in  the  fields  of  environmental  and  clinical  anal- 
ysis. In  many  applications,  traceability  of  quality  control  and  measurement  processes  to 
the  national  measurement  system  is  carried  out  through  the  mechanism  and  use  of  SRM's. 
For  many  of  the  Nation's  scientists  and  technologists,  it  is  therefore  of  more  than  passing 
interest  to  know  the  details  of  the  measurements  made  at  NIST  in  arriving  at  the  certified 
values  of  the  SRM's  produced.  An  NIST  series  of  papers,  of  which  this  publication  is  a 
member,  called  the  NIST  Special  Publication  -  260  Series,  is  reserved  for  this  purpose. 

The  260  Series  is  dedicated  to  the  dissemination  of  information  on  different  phases  of  the 
preparation,  measurement,  certification,  and  use  of  NIST  SRM's.  In  general,  much  more 
detail  will  be  found  in  these  papers  than  is  generally  allowed,  or  desirable,  in  scientific  jour- 
nal articles.  This  enables  the  user  to  assess  the  validity  and  accuracy  of  the  measurement 
processes  employed,  to  judge  the  statistical  analysis,  and  to  learn  details  of  techniques 
and  methods  utilized  for  work  entailing  greatest  care  and  accuracy.  These  papers  also 
should  provide  sufficient  additional  information  not  found  on  the  certificate  so  that  new 
applications  in  diverse  fields  not  foreseen  at  the  time  the  SRM  was  originally  issued  will 
be  sought  and  found. 

Inquiries  concerning  the  technical  content  of  this  paper  should  be  directed  to  the  author(s). 
Other  questions  concerned  with  the  availability,  delivery,  price,  and  so  forth,  will  receive 
prompt  attention  from: 

Office  of  Standard  Reference  Materials 
National  Institute  of  Standards  and  Technology 
Gaithersburg,  MD  20899 

Stanley  D.  Rasberry,  Chief 
Office  of  Standard  Reference  Materials 
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OF  A  SILICON  DIOXIDE  LAYER  ON  SILICON 

G.A.  Candela,  D.  Chandler- Horowitz,  J.F.  Marchiando, 
D.B.  Novotny,  and  B.J.  Belzer 
Semiconductor  Electronics  Division 

M.C.  Croarkin 
Statistical  Engineering  Division 

National  Institute  of  Standards  and  Technology 
(formerly  the  National  Bureau  of  Standards) 
Gaithersburg,  Maryland  20899 

Abstract 

A  Standard  Reference  Material,  SRM-2530,  has  been  designed,  fabricated,  and  certified 
for  the  ellipsometric  parameters  delta,  A,  and  psi,  and  for  the  derived  thickness  and 
refractive  index  of  a  silicon  dioxide  layer  on  silicon  using  a  highly  accurate  ellipsometer 
built  at  NIST.  This  SRM  is  issued  primarily  to  evaluate  the  accuracy  of  ellipsometers.  The 
SRM  consists  of  a  76-mm  (3-in.)  diameter  silicon  wafer  with  a  silicon  dioxide  layer  of  one  of 
three  uniform  thicknesses,  50,  100,  or  200  nm.  The  design  and  fabrication  of  the  SRM  are 
presented  along  with  the  ellipsometric  technique  and  data  analysis  leading  to  certification 
of  this  SRM.  A  least-squares  method  minimizing  the  sum  of  squares  of  deviations  in  A 
and  tf)  between  the  measured  values  and  those  calculated  from  a  model  has  been  used 
in  certifying  the  SRM.  The  derived  values  of  the  thickness  and  refractive  index  may  be 
determined  by  using  either  a  two-layer  or  a  one-layer  model.  The  two-layer  model  assumes 
a  silicon  dioxide  layer  on  a  thin  interlayer  atop  the  silicon  substrate,  whereas  the  one- 
layer  model  assumes  a  single  dielectric  layer  for  the  silicon  dioxide  without  the  interlayer. 
The  two-layer  modeling  analysis  gives  better  agreement  to  the  collective  multiple  sample 
ellipsometric  measurement  data  than  does  the  one-layer  modeling  analysis,  and  gives  a 
value  for  the  refractive  index  of  the  silicon  dioxide  layer  that  is  independent  of  thickness. 
Therefore,  the  certified  values  of  thickness  and  refractive  index  are  based  on  the  two-layer 
model. 


Key  words:  ellipsometric  parameters;  layer  refractive  index;  layer  thickness;  modeling; 
silicon  dioxide  layer;  silicon  substrate;  Standard  Reference  Materials. 
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I.  Introduction 


A  computer-controlled  variable-angle-of-incidence  spectroscopic  ellipsometer  of  high  accu- 
racy has  been  designed  and  constructed  at  NIST.  This  ellipsometer  can  be  used  in  the 
following  modes  of  operation:  conventional  null,  principal  angle  of  incidence,  and  rotating 
analyzer.  The  instrument  is  primarily  used  for  the  metrology  of  semiconductor  materials, 
the  calibration  of  reference  standards,  and  for  research  projects  investigating  multiple- 
layered  structures.  Further  details  regarding  the  ellipsometer  may  be  found  in  reference  1. 

One  of  the  methods  for  transferring  the  capabilities  of  this  instrument  to  the  semiconductor 
industry  is  to  use  this  ellipsometer  to  certify  the  measured  ellipsometric  quantities,  A 
and  -0,  and  the  derived  silicon  dioxide  layer  thickness  and  refractive  index  of  Standard 
Reference  Materials  (SRM's)  [2,3].  The  SRM  can  be  used  as  an  aid  in  the  calibration 
and  the  evaluation  of  optical  and  mechanical  thickness-monitoring  instruments  as  well  as 
ellipsometers  for  which  it  was  specifically  designed  [2,4]. 

The  first  batch  of  SRM-2530's  consists  of  approximately  52  wafers  in  three  thickness  cat- 
egories. The  thickness  categories  of  50,  100,  and  200  nm  refer  to  the  nominal  thicknesses 
of  silicon  dioxide  layers  atop  silicon  substrates.  This  publication  describes  the  design  and 
fabrication  of  the  SRM's  and  the  ellipsometric  and  statistical  methods  for  their  certifica- 
tion. 

Section  II  describes  the  design  and  fabrication  of  the  SRM.  Section  III,  on  ellipsometric 
theory,  defines  the  two  ellipsometrically  measured  quantities,  A  and  ij),  that  are  certified 
at  the  principal  angle  of  incidence  using  the  photometric  method  of  ellipsometry  with  a 
rotating  analyzer.  Sources  of  instrumental  uncertainties  are  assessed  in  section  IV,  and  the 
analysis  of  the  precision  in  the  Fourier  coefficients  that  describe  the  periodically  varying 
light  intensity  is  presented  in  section  V.  Table  I  and  Table  II  summarize  these  results. 
Section  VI  and  Table  III  show  how  the  total  uncertainties  in  A  and  if)  are  derived  from 
systematic  and  random  uncertainties. 

Section  VII  discusses  the  least-squares  procedure  used  to  determine  the  best  values  of  the 
specimen  model  parameters.  In  section  VIII,  the  least-squares  procedure  is  used  to  analyze 
the  two-layered  structure,  consisting  of  a  silicon  dioxide  layer  on  top  of  a  thin  transition 
layer  on  silicon  [5-7].  Section  IX  uses  the  sensitivity  analysis  to  compute  the  uncertainties 
in  the  model  parameters.  Table  IV  shows  the  values  for  the  model  parameters  found  from 
the  measurements  and  their  uncertainties  as  they  would  appear  on  the  certificate  for  a 
typical  SRM  at  each  of  three  different  silicon  dioxide  layer  thicknesses.  In  section  X,  a 
one-layer  model,  where  there  is  no  transition  interlayer  and  each  thickness  group  of  SRM's 
is  allowed  to  have  a  different  layer  refractive  index,  is  presented  and  used  to  compute  the 
corresponding  layer  thicknesses.  In  section  XI,  a  goodness-of-fit  test  is  used  to  compare 
the  two  models.  In  section  XII  and  Appendix  5,  a  computer  program  for  ellipsometric 
calculations  is  described.  This  program  can  be  used  to  calculate  the  layer  thicknesses  for 
the  one-  and  two-layer  models  based  on  certified  A  and  values. 
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II.  Design  and  Fabrication  of  SRM-2530 


Silicon  dioxide  on  silicon  is  one  of  the  better  characterized  layered  structures,  both  chemi- 
cally and  physically.  Thermally  grown  silicon  dioxide  layers  on  silicon  substrates  are  stable 
and  easy  to  clean,  and  extensive  studies  exist  of  the  optical  properties  [8]  for  both  the  layer 
and  bulk  substrate  material.  These  silicon  dioxide  layers  are  the  most  widely  used  and 
measured  layers  in  the  semiconductor  industry.  The  SRM  was  designed  to  help  fulfill  a 
metrological  need  of  the  semiconductor  industry  concerning  the  thickness  and  refractive 
index  of  silicon  dioxide  layers. 


Design  of  SRM 

The  SRM  design  consists  of  a  76-mm  (3-in.)  diameter  silicon  wafer  with  a  silicon  dioxide 
layer  of  uniform  thickness  in  the  shape  of  a  rectangle,  2.5  cm  by  6  cm.  Four  long,  narrow 
(8  mm  by  0.5  mm)  rectangular  windows  were  etched  through  the  dioxide  layer,  and  a 
chromium  overlayer  was  deposited  such  that  its  shape  bisects  the  silicon  dioxide  layer  into 
two  approximately  equal  areas  (see  Fig.  1)  to  be  used  for  the  ellipsometric  measurements. 
The  windows  in  the  silicon  dioxide  are  approximately  half  covered  by  the  chromium  pat- 
tern. The  patterns  of  the  windows  under  the  chromium  can  be  used  for  interferometric 
and  stylus  profilometric  measurements  [2].  The  chromium  overlayer  surface  can  be  useful 
in  aligning  the  ellipsometer  components. 


Silicon  Wafers 

High-quality  silicon  wafers  from  the  same  crystal  were  used  for  the  fabrication  of  each 
batch  of  SRM's.  The  silicon  material  is  boron-doped,  p-type,  (100)  surface,  3  to  10  ohm- 
cm  resistivity.  These  wafers  were  ground  flat  on  the  back  side  and  ground  flat  and  polished 
on  the  front  side.  The  surface  flatness  was  determined  by  using  a  Tropel*  wafer  flatness 
tester.  The  flatness  was  found  to  be  within  one-quarter  wavelength  of  light  at  the  632.8-nm 
helium-neon  laser  line  over  the  measurement  spot  size  of  3-  to  10-mm  diameter  on  the  wafer. 


*  Certain  commercial  equipment,  instruments,  or  materials  are  identified  in  this  paper  to 
adequately  specify  the  experimental  procedure.  Such  identification  does  not  imply  rec- 
ommendation or  endorsement  by  the  National  Institute  of  Standards  and  Technology,  nor 
does  it  imply  that  the  materials  or  equipment  identified  are  necessarily  the  best  available 
for  the  purpose. 
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Fabrication  of  SRM 


The  silicon  dioxide  was  grown  at  a  temperature  of  1050°C  in  a  dry  oxygen  atmosphere. 
The  time  of  growth  was  varied  to  obtain  the  three  thicknesses  of  50,  100,  and  200  nm. 
The  uniformity  of  this  dioxide  layer  thickness  was  measured  using  a  commercial  reflectance 
spectrometer,  and  only  wafers  whose  dioxide  thicknesses  were  uniform,  varying  by  less  than 

0.  6  nm  over  both  of  the  rectangular  areas  of  the  silicon  dioxide,  were  used  for  SRM's.  A 
photolithographic  technique  using  a  photoresist  coating  deposited  on  the  dioxide  and  an 
etchant  of  buffered  hydrofluoric  acid  was  used  to  pattern  the  desired  rectangular  dioxide 
and  window  areas.  The  photoresist  was  then  removed  with  acetone,  and  the  wafers  washed 
with  deionized  water  and  soaked  for  15  minutes  in  hot  concentrated  sulfuric  acid  containing 
freshly  dissolved  ammonium  persulfate.  After  this  treatment,  the  wafers  were  washed  in 
deionized  water  and  spun  dry. 

In  preparation  for  chromium  deposition,  the  wafer  was  immersed  briefly  in  2%  hydroflu- 
oric acid  to  remove  the  native  oxide  from  the  bare  silicon  areas.  A  layer  of  chromium 
approximately  150  nm  thick  was  sputtered  over  the  entire  wafer  surface.  The  chromium 
was  then  patterned  and  etched  with  Cyantek  CR-9*.  The  patterning  resulted  in  leaving 
approximately  half  of  each  of  the  windows  covered  with  chromium  as  is  shown  in  Figure 

1.  All  ellipsometric  measurements  were  made  on  the  two  areas  designated  in  Figure  1  as 
left  side  and  right  side.  The  areas  are  approximately  2.5  cm  by  2.5  cm  located  0.5  cm  to 
the  left  and  right  of  the  chromium  overlayer. 

Cleaning  of  SRM's 

Even  when  the  SRM's  are  stored  in  clean  containers,  a  contaminating  film  forms  on  the 
surface  of  the  silicon  dioxide.  This  unwanted  film  forms  rapidly,  usually  0.1  or  0.2  nm  per 
week  [3],  and  is  removed  before  measurements  are  made  by  rinsing  the  SRM  with  reagent 
grade  ethanol  and  then  with  deionized  water.  Each  solvent  is  blown  off  the  wafer  surface 
with  clean  nitrogen.  The  sample  is  then  permitted  to  stabilize  in  the  atmosphere  for  20  to 
30  minutes  before  measurements  are  made. 


III.  Ellipsometric  Theory 
The  Fresnel  complex  amplitude  reflection  coefficients  are 

H,-f*=|*t«"F.  (1) 


R, 


(2) 


where  Rp  and  Rs  are  the  ratios  of  the  reflected  electric  field  vector  amplitudes,  Erp  or  Ers, 
to  the  incident  electric  field  amplitudes,  Eip  or  i£ta,  for  light  polarized  either  parallel,  p, 
or  perpendicular,  s,  to  the  plane  of  incidence.  The  phase  shifts  of  the  electric  field  induced 
upon  reflection  are  Ap  and  Aa,  for  light  polarized  parallel  and  perpendicular  to  the  plane 
of  incidence,  respectively.  The  ellipsometric  parameters  t/>  and  A  are  defined  by  the  ratio 
of  Rp  to  Rs  [9]: 


^  =  tanVeiA,  (3) 


where 


A=AJ>-A„  (4) 

(5) 


In  the  rotating  analyzer  ellipsometric  (RAE)  method  which  was  used  to  certify  the  SRM's, 
monochromatic  coherent  light  passes  through  a  polarizer  and  onto  the  sample  surface  at 
a  given  angle  of  incidence.  In  this  photometric  method,  the  reflected  beam  goes  through 
a  rotating  analyzer  and  the  light  intensity  is  measured  by  a  silicon  photodiode  detector. 
The  light  intensity,  I,  as  a  function  of  time,  is  [10,11]  (eq.  (Al.l)): 


I(t)  =  J0  (l  +  ot  cos  2A(t)  +  P  sin  2A(t)) ,  (6) 

where  Iq  is  the  average  intensity  and  A(t)  is  the  rotating  analyzer  azimuth.  The  Fourier 
coefficients,  a  and  /?,  are  given  by 


tan2  y>  -  tan2  P 
tan2  V>  +  tan2P 

2  tan  tjj  tan  P  cos  A 
tan2  ij)  +  tan2  P  ' 


ktm  w  —  i an  -i  /_N 


(8) 


where  P  is  the  measured  value  of  the  polarizer  azimuth  with  respect  to  the  plane  of 
incidence. 
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The  Fourier  coefficients,  a  and  /?,  are  obtained  from  a  discrete  Fourier  transform  of  the 
measured  light  intensity  as  a  function  of  time.  Knowing  P  and  the  values  of  a  and  /?,  both 
A  and  if>  may  be  calculated  from  the  following  expressions  ((eqs.  (A1.3),  (Al.4)): 


if)  =  tan  1 


A  =  ±  cos" 


ll  +  a 


tanP 


0  <  fj)  <  90°, 


180°<A  <  180°. 


(9) 
(10) 


The  sign  of  A  must  be  determined  by  a  separate  procedure  such  as  inserting  an  optical 
element  of  known  phase  retardation  between  the  sample  and  rotating  analyzer  [11]. 


IV.  Sources  of  Instrumental  Uncertainties 

This  section  discusses  the  various  uncertainties  arising  from  the  components  of  the  ellip- 
someter.  Estimates  are  made  of  the  uncertainties  in  the  wavelength,  6\]  in  the  angle  of 
incidence,  8<f>;  in  the  polarizer  and  analyzer  azimuths,  6P  and  6 A]  of  effects  due  to  polar- 
izer ellipticity;  and  of  effects  due  to  detector  nonlinearity.  These  systematic  and  random 
components  of  instrumental  uncertainties  are  listed  in  Table  I  and  are  used  in  section  VI 
to  estimate  the  total  uncertainty  in  the  measurement  of  A  and  rf). 

Wavelength  Uncertainty,  6 A 

A  helium-neon  laser  source  is  used  for  both  sample  alignment  and  the  ellipsometric  mea- 
surements. It  is  a  typical  2-mW  laser  whose  linearly  polarized  output  beam  is  converted 
to  circular  polarized  light  and  is  expanded  and  collimated  to  just  over  10  mm.  The  un- 
certainty in  the  wavelength,  632.82  nm,  as  shown  in  Table  I  is  ±0.01  nm  and  is  given  as  a 
systematic  error. 

Goniometer  Calibration  Uncertainty 

The  two  goniometers  that  control  the  azimuth  of  the  polarizer  and  the  angle  of  incidence  of 
the  sample  were  calibrated  to  correct  for  the  difference  between  their  true  angular  positions 
and  their  encoder  readout  dial  positions  [12].  Each  goniometer  is  computer-controlled  by  a 
stepping  motor  and  encoder.  A  calibration  curve  for  each  goniometer  was  determined  and 
stored  in  the  computer  so  that  all  angular  dial  readings  are  corrected  to  within  a  systematic 
component  of  uncertainty  of  ±0.001°.  All  goniometer  angular  positions  are  reached  in  the 
same  direction  of  rotation  to  prevent  an  uncertainty  due  to  the  small  amount  of  backlash 
in  the  goniometer  motor  drive  system. 
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Angle-of-Incidence  Uncertainty,  8<f> 

The  angle  of  incidence  is  determined  by  interferometrically  aligning  the  surface  of  the 
sample  normal  to  the  incident  laser  beam  and  then  using  a  calibrated  goniometer  to  rotate 
the  sample  to  a  given  angle  of  incidence.  This  alignment  procedure  is  as  follows.  First,  the 
surface  of  the  sample  is  aligned  very  close  to  normal  incidence  by  rotating  its  goniometer 
and  recording  the  dial  reading.  Then  fine  adjustments  in  the  two-axis  wafer  stage  are 
made  in  the  horizontal  angular  position,  affecting  the  angle  of  incidence,  and  also  in  the 
vertical  tilt,  affecting  the  polarizer  and  analyzer  azimuths,  to  null  or  minimize  the  number 
of  fringes  viewed  by  means  of  a  beam  splitter  permanently  mounted  between  the  laser 
and  polarizer  element  and  a  TV  monitor.  These  fringes  are  caused  by  the  interference 
between  the  reflected  laser  beam  from  the  wafer's  surface  and  the  reflected  return  beam 
from  the  laser's  exit  window  surface.  This  alignment  at  normal  incidence  is  accurate  and 
reproducible  to  within  ±0.002°  in  the  plane  of  incidence  and  also  ±0.002°  in  the  tilt  of 
the  plane  of  incidence  as  shown  in  Table  I.  The  tilt  uncertainty  causes  an  uncertainty 
in  azimuth  alignment  of  both  polarizer  and  analyzer.  The  angle-of-incidence  uncertainty 
results  in  a  systematic  error  for  a  particular  single  set  of  measurements  on  a  wafer,  but  for 
repeated  mountings  of  a  given  sample  or  for  a  set  of  samples,  it  is  a  random  error.  The 
random  uncertainty  associated  with  the  beam  direction  stability  of  the  laser  was  found  to 
be  neglible  for  this  ellipsometer. 

There  are  three  other  components  of  systematic  uncertainties  in  the  angle  of  incidence.  As 
stated,  the  calibration  of  the  goniometer  produces  a  systematic  uncertainty  of  ±0.001°. 
In  addition,  the  initial  beam  deviation  through  the  polarizer  in  transmission,  0.010°,  was 
corrected  to  be  less  than  ±0.001°  with  the  addition  of  an  optical  wedge  plate  attached  to 
the  polarizer  on  the  entrance  beam  side.  The  uncertainty  associated  with  beam  divergence, 
because  of  collimation,  is  negligible. 

Analyzer  Encoder  Pulse  Uncertainty 

The  analyzer  prism  is  housed  in  a  hollow  shaft  optical  encoder.  This  incremental  encoder 
gives  the  location  of  the  analyzer  azimuth  that  is  rotating  at  approximately  50  rpm.  It  has 
a  once-per-revolution  zero  reference  or  trigger  pulse  in  addition  to  4096  uniformly  spaced 
incremental  pulses  that  are  accurate  to  ±0.001°.  A  divide-by-eight  circuit  is  used  to  allow 
512  pulses  per  revolution  to  trigger  the  capture  of  the  intensity  signal  as  a  function  of  the 
analyzer's  position. 

Azimuth  Alignment  Uncertainties,  6PC  and  6AC 

The  calcite  polarizer  and  analyzer  prism  elements  are  of  the  Glan-Thompson  type  and 
have  an  extinction  coefficient  of  approximately  10~6.  These  polarizing  prisms  must  be 
calibrated  to  determine  the  correction  necessary  to  place  their  transmissive  axes  into  the 
same  zero  azimuth  as  defined  by  the  plane  of  incidence.  This  procedure  is  described  in 
detail  in  references  10  and  11  and  is  outlined  here.  The  true  polarizer  azimuth  position, 
P,  and  analyzer  azimuth  position,  A,  are  obtained  from  their  respective  goniometer  or 
encoder  azimuth  positions  by  adding  to  that  reading  the  correction  or  offset  component, 
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Pc  or  Ac.  The  offset  in  the  analyzer  results  in  a  phase  shift  of  the  intensity  signal  at  the 
detector.  Thus  the  Fourier  coefficients,  a  and  (3  in  eq.  (6),  are  transformed  from  the  Fourier 
coefficients  computed  from  the  measured  intensity  signal.  Calibration  measurements  are 
taken  to  determine  these  offset  angles  for  the  polarizer  transmissive  axis  azimuth  and  its 
goniometer's  zero  dial  reading  and  the  analyzer  transmissive  axis  azimuth  and  its  zero- 
point  optical  encoder  signal.  These  measurements  were  performed  within  a  few  degrees 
of  the  principal  angle  of  incidence  of  each  wafer  to  be  characterized.  Performing  the 
calibration  measurements  at  this  angle  of  incidence  minimizes  the  effect  of  residual  polarizer 
ellipticity  [11].  Each  measurement  procedure  involves  sweeping  the  polarizer  ±2.5°  about 
its  estimated  0°  azimuth  position  in  steps  of  0.5°,  and  at  each  position  obtaining  values 
for  the  measured  Fourier  coefficients.  A  calibration  is  made  by  fitting  the  values  of  each 
expression,  (a2  +  j32)  and  tan-1  (a//?),  versus  the  polarizer  position  where  a  and  (3  are 
the  averaged  coefficients  obtained  from  discrete  Fourier  transforms  on  the  intensity  signal 
(eq.  (A1.12)).  Repetitions  of  the  calibration  measurement  are  then  analyzed  to  determine 
the  mean  values  of  Pc  and  Ac,  and  the  resulting  standard  deviations  are  used  to  evaluate 
the  offset  uncertainties,  6PC  and  6AC.  Typical  values  for  both  6AC  and  8PC  are  ±0.002°. 
These  are  listed  in  Table  I  as  systematic  uncertainties  for  a  single  measurement  but  can 
become  random  uncertainties  for  repeated  measurements. 

Polarizer  Ellipticity 

The  polarizer  ellipticity  was  estimated  to  be  approximately  0.010°  from  many  measure- 
ments made  in  two  zones,  polarizer  azimuth  settings  of  +P  and  —  P  [13].  This  ellipticity 
is  typical  for  a  good  polarizer  (ref.  9,  p.  389).  After  zone  averaging,  the  effects  of  the 
polarizer  ellipticity  become  negligible. 

Detector  Nonlinearity  and  Other  Intensity  Signal  Effects 

The  detector  used  with  the  helium-neon  laser  for  the  intensity  measurements  is  a  sili- 
con photodiode.  The  nonlinear  response  of  this  detector  and  its  amplifying  circuit  [14], 
although  very  small,  was  minimized  by  working  near  the  principal  angle  of  incidence 
and  having  the  polarizer  azimuth  close  to  ^,  in  which  case  the  alternating  component 
of  detected-light  intensity  is  near  zero.  Careful  beam  alignment  with  the  axis  of  rotation 
of  the  analyzer  prism  ensured  a  signal  that  contained  only  the  dc  and  second  harmonic 
signal.  The  background  signal  level,  obtained  when  the  laser  light  intensity  is  extinguished, 
was  found  to  be  negligible  compared  to  the  resolution  of  the  detection  system. 
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Table  I 

Instrumental  Components  of  Uncertainty 


Quantity 

Wavelength  uncertainty,  8X 

Angle-of-incidence  uncertainty,  8(j> 
Sources  of  uncertainty: 

<f>  goniometer  calibration 

Sample  alignment  (in  plane  of  incidence) 

Polarizer  beam  deviation 

Polarizer  azimuth  uncertainty,  8P 

Sources  of  uncertainty: 

P  goniometer  calibration 

Sample  alignment  (tilt  of  plane  of  incidence) 

Transmissive  axis  azimuth  alignment,  8PC 

Analyzer  azimuth  uncertainty,  8  A 
Sources  of  uncertainty: 

Sample  alignment  (tilt  of  plane  of  incidence) 

A  encoder  pulse  uncertainty 

Transmissive  axis  azimuth  alignment,  8AC 

Polarizer  ellipticity  (after  zone  averaging) 

Detector  nonlinearity  (at  principal  angle  of  incidence) 


Systematic  Random 
±0.01  nm 


±0.001° 


±0.001° 


±0.001° 


±0.002° 


±0.001° 
±0.002° 

negligible 

negligible 


±0.002c 


±0.002° 


±0.002e 


V.  Analysis  of  the  Fourier  Coefficients 


The  zone-averaged  Fourier  coefficients,  a'  and  /?' ,  are  defined  by  eqs.  (A1.20)  and  (A1.21) 
in  Appendix  1.  Both  the  short-term  and  long- term  precision  were  evaluated  by  computing 
the  standard  deviations  of  the  zone-averaged  Fourier  coefficients.  The  short-term  repeata- 
bility is  equivalent  to  the  resolution  of  the  instrument,  and  the  long-term  repeatability  is 
the  major  contribution  to  the  uncertainty  in  the  measurement. 

Short-term  precision  and  resolution 

The  precision  or  resolution  of  measurement  for  the  rotating  analyzer  ellipsometric  method 
is  obtained  from  short-term  repeated  measurements  taken  on  an  SRM  sample.  Standard 
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deviations  of  zone-averaged  values,  a'  and  /?',  give  the  limiting  resolution  of  the  instru- 
ment. The  random  fluctuations  of  the  light  source  and  the  detector  response,  which  occur 
during  the  time  it  takes  to  make  a  single  determination,  and  the  resolution  of  the  12-bit 
analog-to-digital  converter  affect  these  standard  deviations.  The  time  duration  for  each 
zone  measurement  consisting  of  32  revolutions  of  the  analyzer  is  approximately  0.5  s,  in 
which  time  the  512  points  per  revolution  times  32  revolutions  analog-to- digital  intensity 
conversions  take  place.  Ten  consecutive  zone-averaged  measurements  of  a'  and  /3'  can 
take  approximately  20  minutes  because  of  the  time  it  takes  the  stepper  motor  to  move  the 
polarizer  from  zone  one  to  zone  two.  The  standard  deviations  of  these  measurements,  8a' 
and  8/3' ,  made  on  an  SRM  sample  are  listed  in  Table  II. 

Long-term  precision  and  uncertainty 

In  the  long-term  measurements,  each  sample  was  mounted,  cleaned,  and  measured  on 
both  the  right  and  left  sides.  This  procedure  was  repeated  an  average  of  5  times  for  each 
sample.  These  long-term  measurements  were  performed  on  6  samples,  2  each  from  the 
three  thickness  groups,  50,  100,  and  200  nm  of  silicon  dioxide.  A  total  of  62  measurements 
were  made  over  a  period  of  1  month.  The  measurements  were  made  at  the  same  angle 
of  incidence  and  polarizer  settings  for  samples  of  each  thickness.  The  variations  in  a' 
and  /?'  may  include  fluctuations  caused  by  slight  variation  in  cleaning  the  surface,  sample 
remounting  angle-of-incidence  fluctuations,  measurement  area  positioning  uncertainties, 
and  possible  instrumental  drift  due  to  warm-up  time.  The  long-term  standard  deviations, 
sai  and  spi,  of  a'  and  j3'  are  used  to  determine  the  uncertainties  in  A  and  if>  for  SRM 
certification. 


Table  II 

Short-  and  Long-Term  Precision  in  the  Fourier  Coefficients 

Radians        Degrees        Radians        Degrees  DF" 
6a!  6/3' 
Short-term  precision        0.0001  0.006°  0.0001  0.006°  9 

Long-term  precision         0.0012  0.069°  0.0006  0.034°  50 


*  degrees  of  freedom 
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VI.  Uncertainties  in  the  Ellipsometric  Parameters 


The  uncertainties  in  the  zone- averaged  quantities  A'  and  ip'  depend  upon  various  con- 
tributions, such  as  uncertainties  in  the  Fourier  coefficients,  8  a'  and  6(3' \  in  the  polarizer 
azimuth,  6P\  in  the  angle  of  incidence,  6(f>;  and  in  the  wavelength,  6X.  The  following  ex- 
pressions for  the  uncertainty  in  A'  and  if)'  due  to  6  a' ,  6/31 ',  and  the  systematic  uncertainty 
in  8P  are  shown  from  eqs.  (Al.24)  and  (Al.25)  in  Appendix  1  to  be 


and 


6 A'  < 


1 

1  -  a'  cos  2P 
1 


\6P  y/l 


8a' 


sin2P 
2v/l  -  a'  2 


y/l  -  a' 


,  a'(3' 

£   

1-a' 


(11) 
(12) 


When  measurements  are  taken  close  to  the  principal  angle  of  incidence  [15],  <^p,  where 
A  =  90°,  and  P  —  if),  the  light  intensity  as  a  function  of  the  analyzer  angle  has  minimal 
fluctuation  and  a'  and  (3'  are  both  near  zero.  At  these  settings,  any  small  nonlinearity  of 
the  detector  does  not  contribute  an  error  to  the  a'  and  /?'  terms.  Also,  because  a'  ~  0  and 
(3'  ~  0,  and  because  |sin2P|  <  1,  eqs.  (11)  and  (12)  reduce  in  first  order  to  the  following: 


6a° 


and 


SA'  <  \6f3'\ 


+  \&P\ 


(13) 


(14) 


The  measurements  are  taken  near  the  principal  angle  because  the  values  of  8tj)'  and  £A' 
are  at  a  minimum. 


All  systematic  and  random  contributions  to  the  total  uncertainties  in  A'  and  -0'  are  calcu- 
lated from  eqs.  (Al.28)  and  (A1.29)  and  are  listed  in  Table  III.  The  effect  of  the  systematic 
component  of  uncertainty  in  the  wavelength  is  negligible.  The  systematic  component  of 
uncertainty  in  the  angle  of  incidence  is  taken  to  be  ±0.002°,  arising  from  the  <f>  goniometer 
calibration  and  polarizer  beam  deviation  contribution  of  uncertainties  as  given  in  Table 
I.  Random  sample  alignment  variations  affect  the  long-term  repeated  measurements  and 
therefore  are  included  only  in  the  standard  deviations  of  a'  and  /?'.  The  uncertainties  in 
A'  and  arising  from  the  angle  of  incidence  uncertainty  are  listed  as  upper-limit  values 
for  all  three  SRM  thickness  categories. 

The  polarizer  azimuth  systematic  component  of  uncertainty,  ±0.003°,  arises  from  both 
the  goniometer  calibration  and  transmissive  axis  alignment  contribution  of  uncertainties 
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as  shown  in  Table  I.  Equations  (13)  and  (14)  give  the  contribution  to  the  uncertainties  6 A' 
and  8ij)'  as  given  in  Table  III.  The  long-term  standard  deviations,  sai  and  s^*,  from  Table 
II  are  used  with  eqs.  (13)  and  (14)  to  obtain  the  contribution  of  the  Fourier  coefficients  to 
the  uncertainties,  listed  in  Table  III. 


Table  III 

Systematic  and  Random  Contributions  to  the  Uncertainties  in  A'  and  if)' 

Systematic        Random           6  A'  61/)' 

Wavelength,  6\                     ±0.01  nm                            ±0.000°  ±0.000° 

Angle  of  incidence,  6<f>           ±0.002°                                ±0.008°  ±0.001° 

Polarizer  azimuth,  6P           ±0.003°    ±0.003° 

sa,  (0.0012  radians)                                     0.069°    ±0.035° 

sp,  (0.0006  radians)                                     0.034°  ±0.034°   

6A'T  *  Wt  * 

Total                                                                           ±0.042°  ±0.039° 


*The  total  uncertainties,  SA't  and  Sij)'T,  for  each  wafer  are  reported  on  the  certificate  that 
accompanies  each  SRM  as  one-standard-deviation  limits  to  random  error  plus  systematic 
error;  the  one-standard-deviation  limits  to  random  error  correspond  to  a  confidence  coef- 
ficient of  0.68.  The  prime  superscripts  have  been  omitted  on  the  certificate. 
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VII.  Description  of  the  Least-Squares  Analysis 


In  order  to  go  from  ellipsometer  parameters  A  and  to  material  parameters  of  thickness 
and  index,  the  data  have  to  be  processed  through  a  model.  The  least-squares  problem 
of  comparing  experimental  measurement  with  a  model  is  formulated  in  Appendix  2.  The 
following  expression,  eqs.  (A2.3)  and  (A2.4), 


5jj?£-[(AJ-Atf+«J-tf),]1  (15) 

involving  the  deviations  between  the  measured  values,  ip'  and  A',  and  the  modeled  values, 
if)"  and  A",  is  minimized.  The  sum  is  over  M  total  measurements  including  different 
samples  and  angles  of  incidence.  This  least-squares  analysis  was  applied  to  the  two-layer 
and  one-layer  models. 

VIII.  Parameters  of  the  Two-Layer  Model 

The  following  least-squares-fitting  procedure  was  used  in  determining  the  parameters  for 
the  SRM's.  In  analyzing  the  two-layer  model,  all  the  data  obtained  on  all  the  samples 
for  the  three  silicon  dioxide  thicknesses,  50,  100,  and  200  nm,  were  analyzed  collectively. 
As  stated  in  section  II,  all  the  wafers  were  obtained  from  the  same  silicon  ingot,  and 
all  the  silicon  dioxide  films  were  grown  under  the  same  conditions  except  for  the  length 
of  time  for  dry  silicon  dioxide  growth.  The  layer  parameters  that  are  varied  in  the  least- 
squares  analysis  are  determined  when  |g|  is  at  a  minimum.  At  this  minimum,  the  deviations 
between  the  measured  and  modeled  A's  and  should  be  on  the  order  of  the  measurement 
uncertainties  in  A'  and  tp' . 

A  depiction  of  the  two-layer  model  is  shown  in  Figure  2.  The  model  involves  a  thin  in- 
terlayer  between  the  silicon  substrate  and  silicon  dioxide  layer.  Two  parameters,  common 
to  all  samples,  that  were  not  estimated  from  the  measurement  data  were  the  imaginary 
part  of  the  silicon  substrate  refractive  index,  ka  =  0.018  [8],  and  the  interlayer  refractive 
index,  n;  =  2.8  [5-7].  The  model  parameters  that  were  estimated  by  least  squares  include 
the  silicon  dioxide  layer  thickness,  for  each  sample  and  the  following  parameters,  com- 
mon among  all  samples:  the  silicon  dioxide  layer  refractive  index,  n/;  the  interlayer  film 
thickness, -U;  and  the  real  part  of  the  silicon  substrate  refractive  index,  n,.  The  imaginary 
part  of  the  silicon  dioxide  layer  refractive  index  at  632.8  nm  is  zero,  and  it  was  found  that 
setting  the  imaginary  part  of  the  interlayer  refractive  index  to  zero  had  negligible  effect  on 
the  least-squares  fit. 

From  a  total  of  239  ellipsometric  measurements  on  52  samples,  including  left  and  right  side 
and  long-term  repeated  measurements,  the  best  estimates  from  the  least-squares  analysis 
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Table  IV 

Parameters  and  Uncertainties  for  the  Two-Layer  Model* 


Values 

from  the  Least- Squares  Analysis 

Sample  Serial  No. 
390-03         384-48  389-10 
tf               tf  tf 

u 

All  Samples 

nf  n„ 

Igl  a>b>c 

=  0.068° 

53.9 

97.9 

202.3 

1.0 

1.461 

3.875 

Contributions  to  Total  Uncertainties  from  the 

Sensitivity  Analysis 

6tf 

6tf 

6tf 

6U 

6nf 

6n, 

6k,  f  = 
6m  f  = 

0.078° 
0.002° 

±0.005 

±0.4 

±0.11 
±0.02 
±0.01 
±0.06 

±0.13 
±0.01 
±0.01 
±0.06 

±0.18 
±0.01 
±0.01 
±0.06 

±0.03 
±0.01 
±0.09 
±0.21 

±0.0003 
±0.0000 
±0.0000 
±0.0000 

±0.0011 
±0.0001 
±0.0001 
±0.0010 

Total  « 

±0.18 

±0.20 

±0.25 

±0.33 

±0.0003 

±0.0022 

Final  Values  and  Their  Uncertainties* 

tf  +  u 

tf+U 

tf  +  U 

U 

nf 

n. 

54.9 

98.9 

203.3 

1.0 

1.461 

3.875 

6{if+U) 

6(tf+U) 

6(tf+U) 

su 

6nf 

6na 

±0.5 

±0.5 

±0.6 

±0.3 

±0.001 

±0.003 

*  Thicknesses  in  nanometers 

aSee  eq.  (15) 

hka  =  0.018,  see  ref.  8 

crij  =  2.8,  see  refs.  5,  6,  and  7 

d  Random  component  of  uncertainty,  see  eq.  (A3. 8) 

e  Random  component  of  uncertainty  (Systematic  component  of  ±0.002°  gives  negligible 
effect) 

'  Systematic  components  of  uncertainty 

9  See  eq.  (A3. 7)  for  the  calculation  of  the  total  uncertainties 

*The  total  uncertainties  for  each  wafer  are  reported  on  the  certificate  that  accompanies 
each  SRM.  The  uncertainties  are  computed  as  one-standard-deviation  limits  to  random 
error  plus  systematic  error;  the  one-standard-deviation  limits  to  random  error  correspond 
to  a  confidence  coefficient  of  0.68. 


14 


for  the  model  parameters  are  shown  in  the  first  row  of  numbers  in  Table  IV.  Three  repre- 
sentative samples  (Sample  Nos.  390-03,  384-48,  and  389-10)  are  shown.  The  value  for  the 
least-squares  parameter,  |g|,  is  0.068°.  The  following  common  parameter  values  are  found: 
interlayer  thickness  t{,  1  nm;  the  refractive  index  of  the  silicon  dioxide  film  nf,  1.461;  and 
the  value  of  the  silicon  substrate  refractive  index  n4,  3.875. 

IX.  Sensitivity  Analysis  and  Uncertainties  in  the  Model  Parameters 

In  Appendix  3,  a  sensitivity  analysis  of  the  model  parameters  is  discussed.  The  manner  in 
which  one  may  associate  an  uncertainty  with  the  model  parameters  is  developed  leading 
to  eq.  (A3. 7).  Using  the  sensitivity  analysis,  the  value  of  sg  =  0.078°,  the  standard 
deviation  in  the  angle  of  incidence,  3$  —  0.002°  (random  part  of  uncertainty  in  Table  I), 
the  systematic  part  of  the  angle  of  incidence  uncertainty,  6<f)tya  =  0.002°,  an  uncertainty 
of  6k8  =  ±0.005  [8,16]  for  the  value  of  k8  —  0.018,  and  the  uncertainty  in  the  interlayer 
refractive  index  of  Srii  =  ±0.4  [5-7]  for  the  value  of  m  —  2.8,  the  relative  contributions  to 
the  total  uncertainties  in  the  model  parameters  can  be  found. 

The  second  row  of  Table  IV,  labeled  a5,  gives  random  components  of  uncertainties  in  the 
estimated  model  parameters  due  to  the  fit  of  the  measured  data.  The  next  row  labeled 
gives  the  contribution  of  uncertainty  due  to  the  random  component  of  angle  of  incidence 
uncertainty  (systematic  component  of  ±0.002°  gives  negligible  effect).  Rows  labeled  8ka 
and  6ni  give  the  respective  contribution  to  the  systematic  components  of  uncertainties. 

The  value  of  this  sensitivity  analysis  is  that  the  contributions  to  the  uncertainties  can  be 
pinpointed.  The  two  largest  sources  of  uncertainty  in  the  thickness  of  the  total  oxide  arise 
from  the  value  of  sg  obtained  from  the  measurement  fit  and  from  the  large  uncertainty  in 
the  thickness  of  the  interlayer  due  to  Srii.  The  large  uncertainty  in  the  interlayer  thickness 
arises  from  the  correlation  of  the  refractive  index  and  the  thickness  for  this  very  thin  film. 

X.  Parameters  of  the  One- Layer  Model 

Two  different  one-layer  models  have  been  analyzed.  The  first  model  assumes  a  common 
value  of  refractive  index,  independent  of  layer  thickness,  to  all  SRM  samples.  The  second 
model  allows  a  different  common  value  of  refractive  index  for  each  layer  thickness  group 
of  SRM's.  The  experimental  fit  to  the  first  model  is  very  poor  (see  sec.  XI).  The  second 
model,  which  gives  a  layer  refractive  index  that  is  different  for  each  of  the  three  thickness 
categories  of  SRM's,  50,  100,  and  200  nm,  is  not  physically  realistic,  but  is  commonly  used 
when  processing  single-sample  data  and  is  therefore  given  as  supplementary  information 
on  the  certificate. 

The  three  thickness  categories  of  SRM  wafers,  50,  100,  and  200  nm,  were  analyzed  sep- 
arately by  the  least-squares  procedure  using  this  second  single-layer  model.  All  SRM's 
of  the  same  thickness  group  were  allowed  to  have  a  best-fit  common  refractive  index,  n/, 
of  the  film,  and  a  best-fit  common  substrate  refractive  index,  n#,  and  a  fixed  common 
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value  for  k„.  An  individual  thickness  was  calculated  for  each  wafer.  The  value  of  na  was 
found  to  be  the  same  as  that  of  the  two-layer  model,  ns  =  3.875,  and  the  corresponding 
film  refractive  indices,  n/  for  the  50-,  100-,  and  200-nm  wafers,  were  1.468,  1.465,  and 
1.458,  respectively.  For  the  samples  in  Table  IV,  the  one-layer  model  thicknesses  for  serial 
numbers  390-03,  384-48,  and  389-10  are  54.4,  98.1,  and  204.0  nm,  respectively.  These  film 
thicknesses  should  be  compared  with  the  two-layer  calculation  that  includes  the  thickness 
of  the  interlayer  as  given  under  "Final  Values  and  Their  Uncertainties"  in  Table  IV.  The 
two-layer  model  and  the  one-layer  model  for  the  50-,  100- ,  and  200-nm  samples  disagree 
by  the  amounts  0.5,  0.8,  and  0.7  nm,  respectively. 

XI.  Goodness-of-Fit  Test 

Three  distinct  models  were  considered  and  discussed  earlier  in  the  text,  i.e.,  two  one-layer 
models  and  a  two-layer  model.  The  one-layer  model  considered  in  this  section  assumes 
a  common  value  of  layer  refractive  index,  independent  of  thickness.  The  two-layer  model 
has  an  interlayer  between  the  top  silicon  dioxide  film  and  silicon  substrate.  The  model 
parameters  of  this  interlayer  include  its  refractive  index  (nj)  and  thickness  (tfj).  Hence,  the 
two-layer  model  requires  two  parameters  more  than  the  one-layer  model. 

It  is  important  to  test  the  fit  of  the  experimentally  measured  ellipsometric  parameters  A' 
and  tp'  to  those  predicted  by  the  nonlinear  model  used  for  estimating  individual  specimen 
thicknesses  and  indices  of  refraction.  This  is  done  independently  for  each  model  using  a 
goodness-of-fit  procedure  that  is  outlined  in  Appendix  4.  An  F  statistic  for  testing  for 
lack  of  fit  is  given  by  eq.  (A4.4), 

jp         [SStotal  error      SSpvre  error]  /  [df total  error      dfpure  error]  (16) 

[S SpUre  error  ]  I  [dfpure  error  I 

where  the  SS  stands  for  the  sum  of  squares  of  the  errors  (see  eqs.  (A4.1)  and  (A4.3))  and 
df  are  the  degrees  of  freedom  as  defined  in  Appendix  4.  The  F  statistic  is  compared  with 
Fa.{y\,V2)  which  represents  the  upper  a  percent  point  of  the  F  distribution  with  v\  and 
V2  degrees  of  freedom.  The  F  test  is  not  exact  because  the  underlying  model  is  nonlinear; 
however,  a  value  of 

F<Fa{uuu2),  (17) 

indicates  that  there  is  no  significant  lack  of  fit  to  the  collective  data  for  the  model  being 
considered. 

The  pure  error  term  is  constructed  from  repetitions  on  six  specimens  that  were  remeasured 
several  times  at  the  same  wavelength  and  the  same  angle  of  incidence  over  the  course  of 
the  calibration  of  the  batch  of  SRM's.  The  degrees  of  freedom  associated  with  the  pure 


error  term  is  =  92.  The  total  error  term  is  constructed  from  the  deviations  from  the 
fit  with  377  degrees  of  freedom  for  the  two-layer  model;  the  sum  of  squares  for  lack  of  fit 
has  degrees  of  freedom  equal  to  the  degrees  of  freedom  for  total  error  minus  the  degrees  of 
freedom  for  pure  error,  or  v>i  =  285. 

The  value  of  F  from  eq.  (16)  is  3.8  for  the  two-layer  model.  For  a  significance  level  of 
a  =  0.01,  the  corresponding  critical  value  for  .F.oi(285, 92)  is  1.5.  The  F  statistic,  which  is 
very  sensitive  to  lack  of  fit  or  large  outliers,  exceeds  the  critical  value  but  not  by  enough  to 
indicate  an  egregious  lack  of  fit.  The  goodness-of-fit  test  indicates  that  the  two-layer  model 
adequately  represents  the  wafer's  surface.  The  F  statistic  is  21.4  for  the  corresponding 
one-layer  model  or  five  times  greater  than  for  the  two-layer  model,  indicating  a  much 
greater  degree  of  lack  of  fit  to  the  model. 

XII.  Computer  Program  for  Determining  the  Thickness  of  the  SRM  Silicon  Dioxide  Layer 

Accompanying  this  SRM-2530  is  a  floppy  disk  formatted  on  an  IBM  PC/XT  that  contains 
a  computer  program.  This  program  is  a  combination  of  McCrackin's  method  [17]  followed 
by  a  least-squares  method  for  calculating  the  thickness  of  the  SRM  silicon  dioxide  layer.  In 
Appendix  4,  a  description  is  included  of  how  to  use  this  program  to  calculate  the  thickness 
for  both  the  one-  and  two-layer  models. 
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Silicon  Dioxide      tf,  nf 
Layer 

=    Interlayer  tj,  nv  kj 

/     Silicon  ns,  ks 

Substrate 

/ 

Figure  2.  Schematic  of  two-layer  model.  Cross-section  schematic  of  silicon  dioxide 
on  silicon  designating  seven  of  the  parameters  for  the  two-layer  model.  The  one- 
layer  model  assumes  that  no  interlayer  is  present,  t{  =  0. 
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Appendix  1 
Ellipsometer  Precision  and  Uncertainty- 
All  ellipsometric  measurements  on  the  Standard  Reference  Materials  were  obtained  by 
an  ellipsometer  using  the  rotating  analyzer  configuration.  Here,  a  simple  discussion  is 
presented  regarding  the  precision  and  uncertainty  of  the  instrumentation  in  order  to  ob- 
tain the  uncertainty  in  the  two  ellipsometric  parameters,  A  and  ip.  An  assessment  is 
made  of  the  various  random  and  systematic  uncertainties  which  contribute  to  the  total 
uncertainty  of  A  and  ij).  In  this  appendix,  no  assessment  is  made  regarding  the  effects 
due  to  sample  preparation,  mounting,  or  alignment,  e.g.,  relocating  the  illuminated  spot 
on  the  sample  during  remounting  or  those  effects  due  to  cleaning  the  surface  of  the  sam- 
ple. These  issues  are  discussed  in  section  IV. 

For  the  case  of  the  rotating  analyzer  ellipsometer,  the  intensity  flux  of  light,  I(t),  at  the 
detector  is  a  known  function  of  time  and  is  given  by 

I(t)  =  J0(l  +  acos2A(t)  +/?sin2A(<)),  (.41.1) 

where  A(t)  refers  to  the  azimuth  angle  of  the  transmissive  axis  of  the  analyzer,  is  mea- 
sured with  respect  to  the  plane  of  incidence,  is  a  function  of  time,  <,  and  is  given  by 

A(t)  =  2tt  ft  +  Ac. 

Here,  IQ  is  the  average  intensity,  a  and  /3  are  the  normalized  Fourier  coefficients,  /  is 
the  mechanical  rotation  frequency,  and  Ac  is  a  constant  phase  factor.  The  two  Fourier 
coefficients  represent  contributions  of  the  second  harmonic  to  the  measured  intensity. 
Also,  the  measured  intensity  is  known  to  be  a  function  of  the  polarizer  azimuth,  P,  so 
that  the  phase  and  amplitude  information  of  the  reflected  light,  i}>  and  A,  may  be  ex- 
pressed by  [Al.l] 

0  <  -0  <  90°  (.41.3) 


180°  <  A  <  180°  .  (41.4) 


U1.2 


V>  =  tan  1 


and 


A  =  ±  cos  1 


[1  +  ^ 


tanP 
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In  order  to  determine  a  and  j3  from  experiment,  the  intensity  flux  is  sampled  at  512  dis- 
crete equally  spaced  points  of  azimuth  angle  per  mechanical  revolution  of  the  analyzer. 
The  analyzer's  relative  azimuth  angle  is  obtained  from  encoder  pulses  triggered  by  a  zero 
reference  signal  from  its  hollow  shaft  optical  encoder.  Data  are  collected  for  32  revolu- 
tions of  the  analyzer. 

Within  each  mechanical  revolution  of  the  analyzer,  the  measured  intensities  are  assumed 
to  be  of  the  form 


Ij  —  <lq  +  0-2  cos2^4j  +  b2  sin  2  A,  +  ejk  (A1.5) 

=  a0(l  +acos2A,-  +  /?  sin  2Aj)  +  ejk,  (A1.6) 

where 

Aj  =  2*^=A  i  =  l,2,...,J,  (A1.7) 

where  J  is  equal  to  512  and  K  is  equal  to  32,  and  where  j  indexes  the  discrete  positions 
of  azimuth  angle  of  the  analyzer,  k  indexes  the  revolutions  of  the  analyzer,  and  ejk  is 
due  to  small  fluctuations  of  random  noise. 

A  method  of  least  squares  is  then  used  to  provide  a  best  estimate  of  ao ,  ,  and  b2  •  Let- 
ting Ijk  denote  the  measured  intensity  at  the  jth  value  of  azimuth  angle  and  A;th  revolu- 
tion of  the  analyzer,  a  discrete  Fourier  transform  is  performed  where 

K 


a>ok  =  7  y^-fyb,  ao  =  -j£^2aok,  (A1.8) 

j=l  Jb=l 

2    3  1  K 

a2k  =  j^2T3k  cos2Aj,  a2  =  —  ]T  a2k,  (A1.9) 

3=1  k=i 

1    3  1  K 

b2k  =  j       Ji*  sin2Ai>  ^2  =  ^^  hk,  (ALIO) 


j=l  k=l 
U2k  n  &2Jb 


fit  =  — ,  (^1-11) 
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(A1.12) 


The  standard  deviations  of  a  and  /?  are  estimated  by 


^rE(^-a)2  (-41.13) 

k  =  l 

and 

*=i 

The  statistical  uncertainty  associated  with  determining  the  normalized  Fourier  coeffi- 
cients, d  and  /?,  is  assumed  or  defined  to  be 

8a  —  sa  and  8/3  =  sp.  (A1.15) 

There  is  an  inherent  angular  or  phase  offset  between  the  intensity  detected  at  the  first 
encoded  pulse  and  that  associated  with  the  true  zero  azimuth  of  the  analyzer.  This  is 
due  to  two  distinct  conditions,  one  being  a  mechanical  offset  in  the  mounting  of  the  ana- 
lyzer in  the  optical  encoder,  and  the  other  being  a  phase  lag  in  the  signal  processing  cir- 
cuits. Together,  their  net  effect  is  that  of  inducing  a  simple  coordinate  rotation  onto  the 
analyzer  though  an  effective  angle,  Ac,  and  an  amplitude  reduction,  which  for  this  ellip- 
someter,  was  negligible.  By  following  a  simple  calibration  procedure  [Al.2],  this  angular 
offset  (Ac)  is  determined.  Then,  the  calibrated  normalized  Fourier  coefficients  become 


ac  =  acos2Ac  +  (3sin2Ac,  (A1.16) 

(3C  =  dsin2Ac  +/?cos2Ac,  (41.17) 

and  their  corresponding  uncertainties 

8ac  <  \8acos2Ac\  +  |£/?sin2i4c|  +  2\8Ac(asin2Ac  -  /?cos2Ac)|,  (41.18) 

8flc  <  |£dsin2Ac|  +  |£/?cos2Ac|  +  2|^Ac(d  cos  2AC  -  /3  sin2Ac)|.  (A1.19) 


The  last  term  in  each  of  the  above  two  equations  depends  on  the  uncertainty  in  the  an- 
alyzer offset  value  (8AC),  but  is  negligible  when  compared  to  the  first  and  second  terms 
when  a  and  j3  are  small  and  of  the  order  of  8a  or  6j3  at  the  principal  angle  of  incidence. 
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Furthermore,  these  values  (ac,/?c)  are  determined  at  a  specific  position  setting  (P)  of  the 
polarizer,  and  hence  depend  upon  P.  Also,  it  is  known  that  most  polarizers  are  imper- 
fect and  contain  some  inherent  ellipticity.  So,  to  correct  for  this  systematic  error  as  well 
as  for  azimuth  calibration  uncertainties  in  the  polarizer,  one  performs  a  procedure  called 
zone  averaging.  This  is  accomplished  by  determining  the  Fourier  coefficients  (ac,(3c)  at 
the  positive  polarizer  setting  (+P),  which  may  be  denoted  by  (a+,/?+),  as  well  as  at  the 
negative  polarizer  setting  (— P),  denoted  by  (a~ ,/?")•  The  zone-averaged  values  are  de- 
fined by 

i  («++«") 


and 


(/?+-/n 


(41.20) 


(41.21] 


The  negative  sign  present  in  eq.  (Al.21)  is  due  to  the  fact  that  (3  is  an  odd  function  of 
the  polarizer  azimuth  P,  as  may  be  seen  from  eq.  (8).  The  corresponding  uncertainties 


6a'  < 


\6a+\  +  \6a 


and 


60 


t<\^\  +  \^-\m 


(41.22) 


(41.23) 


Knowing  a',  /?',  P,  and  eqs.  (A1.3)  and  (A1.4)  ,  one  may  then  calculate  zone-averaged 
estimates,  A'  and       of  A  and  ijj.  Accordingly,  it  also  follows  that  the  uncertainties  may 

be  estimated  by 


and 


;'  < 


6 A'  < 


1-  a'  cos  2P 
1 

y/l  -  a'  2  -  (3' 


6P  yf\  -a'  2 
,  a'/?' 


6a' 


sin2P 


2y/l  -  a'  2 


6a' 


I -a' 


(41.24) 


(41.25) 


For  ellipsometric  data  taken  at  the  principal  angle  of  incidence,  A'  =  ±90°,  and  with 
the  polarizer  azimuth,  P  =  if)' ,  a'  and  /3'  are  small  and  eqs.  (A1.24)  and  (A1.25)  reduce 

to 


6i/>'  <  ^|sin2P|  +  6P  < 


6a' 


(41.26) 


and 


24 


6 A'  <  \6(3'\. 


(A1.27) 


As  may  be  inferred  from  the  above  equations,  it  is  often  the  case  that  the  uncertainties, 
6  A'  and  6if)' ,  achieve  a  minimum  whenever  measurements  are  taken  near  the  principal 
angle  of  incidence. 

The  values  of  A'  and  if}'  are  determined  simultaneously  at  a  given  angle  of  incidence, 
and  wavelength,  A.  Therefore,  some  accounting  of  uncertainty  in  <f>  and  A  should  be  in- 
cluded within  the  total  uncertainty  in  A'  and  if)' ,  8A'T  and  8if)'T.  Accordingly,  one  may 
formally  write 


and 


tyT  <  w  + 

8j>' 
6(f) 

6(f)  + 

6^' 
6\ 

6A'T  <  6 A'  + 

6A' 

6(f> 

6(f)  + 

6A' 
6X 

sx 


6X. 


(A1.28) 
(A1.29) 


The  incremental  rate  changes  in  A'  and  i'  with  respect  to  (f>  are  determined  experimen- 
tally. These  are  then  multiplied  by  the  uncertainty  6(f)  to  obtain  their  relative  contri- 
bution to  the  total  uncertainty  in  A'  and  i'.  Because  the  magnitude  oi  6X  from  the 
helium-neon  laser  is  particularly  small,  the  incremental  changes  in  A'  and  ip'  with  re- 
spect to  A  may  be  considered  negligible. 
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Appendix  2 


Formulation  of  the  Least-Squares  Problem 


In  order  to  characterize  the  layered  structure  of  the  sample,  it  is  necessary  to  invert  the 
standard  ellipsometric  equations.  The  equations  describe  how  the  material  properties 
of  the  sample  are  able  to  induce  a  phase-shift  in  the  reflected  light  which  is  then  mea- 
sured by  the  ellipsometer.  Because  of  the  complexity  of  the  representation,  it  is  usually 
not  possible  to  find  a  simple  analytic  expression  which  will  invert  the  equations.  One 
common  approach  to  performing  such  an  inversion  is  to  formulate  it  as  a  nonlinear  least- 
squares  problem  [A2.1-4].  Here,  one  considers  a  sequence  of  forward  scattering  problems, 
where  each  increment  of  the  sequence  involves  three  distinct  steps,  where  the  steps  in- 
clude:   starting  with  or  providing  a  good  estimate  of  values  to  the  model  parameters, 
determining  the  deviations  between  the  experiment  and  model,  and  then  updating  the 
model  parameters  with  better  values.  This  sequence  is  continued  or  iterated  until  the 
magnitude  of  the  corrections  becomes  sufficiently  small. 

The  ellipsometric  equation  are  of  the  form: 

A  =  A(<£,A,b)  (42.1) 
V>=y^,A,b),  (42.2) 

where  A  is  the  wavelength  of  the  incident  light,  <j>  is  the  angle  of  incidence,  and  b  is  a 
vector  whose  components  specify  the  model  parameters,  e.g.,  indices  of  refraction,  ex- 
tinction coefficients,  and  thicknesses  for  the  thin  films  and  substrate.  The  standard  pro- 
cedure involves  minimizing  some  non-negative  scalar  error  expression  containing  the 
deviations  between  experimental  and  modeled  values  of  A  and  V>  m  the  least-squares 
sense,  e.g., 

M 

G(b)  =  m  £  N  "  a;,)2  +  W  -  *<)2}  {A2-3) 

i=l 

=  gTg  =  |g|2, 

where  the  prime  denotes  the  value  measured  by  experiment,  the  double  prime  denotes 
the  least-squares  estimate,  M  denotes  the  total  quantity  of  measurements  of  (A,-0)  from 
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experiment,  g  is  a  column  vector  of  the  deviations,  i.e., 

92i-i  =  (A;  -  A?)  /V2M,  (A2.5) 
^  =  (^-^')7\^M,  (A2.6) 
where  1  <  i  <  M,  and  the  superscript  T  denotes  transposition. 

The  third  step  of  the  sequence  is  concerned  with  the  procedure  for  correcting  or  improv- 
ing the  model  parameters  with  better  values,  i.e.,  the  Newton  step.  An  estimation  of  the 
necessary  correction  may  be  realized  by  linearizing  the  functional  representation  of  the 
model  and  solving  the  resulting  matrix  equation,  g=Jv,  where  v  is  a  vector  (Newton 
step)  for  improving  the  model  parameters  that  were  selected  to  undergo  variation,  e.g., 
Vj  ~  6b j,  and  J  is  the  sparse  matrix  involving  the  Jacobian  which  is  not  necessarily 
square,  e.g.,  J{j  ~  (1/v^M)  (dAi/dbj). 

Solving  this  matrix  equation  is  common  to  optimization  problems  involves  only  linear 
algebra,  and  is  well-studied.  It  is  known  that  additional  numerical  stability  or  robustness 
may  result  by  requiring  the  norm  of  v  be  minimized  as  well.  This  can  be  accomplished 
simply  by  modifying  the  error  expression  to 

G  =  (g  -  J  v)T(g  -  J  v)  +  *vT  v,  {A2.7) 

where  k  is  a  positive  scalar  parameter  subjectively  chosen  to  equal  0.01  for  these  calcu- 
lations. It  is  also  known  that  the  columnar  scaling  of  J  affects  the  accuracy  of  the  solu- 
tion and  the  effectiveness  of  k  in  moderating  the  rate  of  convergence  for  finding  the  solu- 
tion v.  Here,  a  simple  choice  for  the  scaling  can  be  found  by  using  the  diagonal  elements 
from  JT  J  and  then  defining  the  diagonal  matrix,  S,  where 

s»  =  v/(jTj);> '  {A2'8) 

Then  letting 

J'=JS-\  (42.9) 

v'  =  Sv,  (A2.10) 

r  =  g-Jv  =  g-jV,  (42.11) 
a  suitable  error  expression  may  be  defined  by 

G  =  (g  -  J'v')T(g  -  J  V)  +  «v'Tv'  ,  (A2.12) 

(A2.13) 
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The  criteria  for  critical  points  or  relative  minima  is  a  vanishing  variation,  i.e., 
dG/dv'j  =  0,  which  yields  a  set  of  equations  that  may  be  expressed  as 


1 

J' 

r 

g 

.J,T 

-k1 

v' 

0 

(42.14) 


that  must  be  solved  for  v'. 

Because  of  the  sparsity  of  the  Jacobian,  it  is  expedient  to  utilize  an  iterative  method  for 
solving  the  above  matrix  algebra  problem.  Algorithms  [A2.5]  exist  that  specifically  ad- 
dress this  type  of  problem,  among  which  one  [A2.6,7]  utilizes  a  relatively  stable  Lanczos 
process  (Krylov  space  decomposition)  in  formulating  the  method  of  steepest-descents. 
Essentially,  the  method  requires  that  each  updating  vector  be  orthogonal  to  the  previous 
update  vectors.  Further  details  may  be  found  elsewhere  [A2.6-10].  From  this,  one  finds 

v',  which  then  leads  to  v  (the  Newton  step),  which  is  then  used  to  improve  the  estimate 

e  i        i       e     li  i_  .new  old 

of  the  values  for  the  model  parameters,  b,  e.g.,  bj       —  b-     +  Vj. 

Using  this  improved  b,  the  sequence  is  iterated  again  until  either:    |  g  |  becomes  suffi- 
ciently small  of  the  order  of  a  few  milli-degrees,  i.e.,  the  resolution  of  the  measurement, 
or  until  |v'|  becomes  sufficiently  small  so  that  |  g  |  suffers  no  further  reduction  regardless 
of  magnitude.  It  is  especially  in  this  last  case  that  it  becomes  necessary  to  scan  a  grid  of 
model  parameters.  Multiple  pseudo- minima  may  be  encountered,  e.g.,  non-uniqueness. 
Often,  this  serves  to  reveal  either:    (i)  correlation  which  prevents  model  parameters  from 
being  resolved  independently,  i.e.,  the  measurement  data  were  not  sufficiently  function- 
ally independent  which  thereby  induces  a  functional  dependence  among  the  model  pa- 
rameters; or  (ii)  the  inadequacy  of  the  model  in  providing  a  sufficiently  good  physical 
description  of  the  process,  which  is  likely  whenever  |  g  |  exceeds  much  beyond  the  resolu- 
tion of  the  measurements. 

Finally,  it  is  important  to  realize  that  in  the  above  outlined  steps,  the  main  empha- 
sis was  that  of  searching  for  and  ascribing  good  values  to  the  parameters  of  a  physical 
model  which  had  already  been  specified.  It  is  possible  that  two  distinct  models  may 
reduce  the  deviations  in  |  g  |  to  an  equivalent  magnitude.  The  problem  of  characteriz- 
ing the  sample  then  becomes  one  of  comparing  models  and  thereby  selecting  the  better 
model.  One  heuristic  approach  has  been  to  select  that  model  which  provides  the  small- 
est deviations  in  |  g  |  while  using  the  fewest  modeling  parameters  and  being  consistent 
with  physical  reality.  But  a  comparison  assumes  an  ordering,  and  that  requires  a  num- 
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ber.  So  the  problem  becomes  one  of  reducing  a  model  to  a  number.  Although  this  re- 
duction is  certainly  nontrivial,  it  is  expedient  to  follow  simple  statistics  and  to  consider 
the  use  of  the  F  statistic  in  assessing  a  so-called  goodness-of-fit  test.  For  further  discus- 
sion on  this  topic,  the  reader  is  referred  to  Appendix  4. 
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Appendix  3 


Sensitivity  Analysis  of  the  Model  Parameters 


During  the  calculation  of  the  least-squares  problem,  the  main  steps  center  on  finding 
accurate  numerical  values  for  the  model  parameters.  It  then  follows  naturally  that  the 
next  logical  step  would  involve  some  assessment  of  uncertainties  associated  with  knowing 
those  values.  Fortunately  for  a  restricted  class  of  linearizable  problems,  one  may  provide 
estimates  of  the  uncertainties  by  utilizing  a  formalism  similar  to  that  used  in  solving  the 
least-squares  problem. 

Following  the  notation  presented  in  Appendix  2,  let  the  functional  representation  of  the 
model  be  expanded  to  first  order  about  the  assumed  fixed  point,  e.g., 


A' -A" 


j=i  3 


N 

E 


Sbj  + 


Sbj  + 


0A, 


drffi 


6<f>, 


S<f>, 


(A3.1) 


(A3.2) 


where  N  denotes  the  quantity  of  distinct  model  parameters,  and  <f>  refers  to  the  angle  of 
incidence.  The  values  (bj)  and  uncertainties  (Sbj)  of  the  model  parameters  were  either 
determined  from  analyzing  the  measurements  or  taken  from  the  literature.  The  litera- 
ture values  and  associated  uncertainties  that  are  assumed  known  and  unchanged  must 
be  considered  when  estimating  the  uncertainty  in  the  model  parameters  estimated  from 
the  measurement  data.  This  leads  to  a  partitioning  of  the  uncertainties  Sbj  into  those 
unchanged  (u)  and  those  varying  (v),  so  that  the  above  expansion  may  be  written  as 


g  =  J„  v+Juu+J^*, 


(A3. 3) 


where  J  denotes  the  Jacobian  appropriate  to  the  partitioning,  and  $  is  a  columnar  array 
involving  the  uncertainty  in  the  angle  of  incidence.  Then,  the  problem  of  determin- 

ing uncertainties  in  the  values  of  the  variant  model  parameters  is  that  of  estimating  the 
least  upper  bound  to  \vj\.  A  formal  solution  of  the  above  equation  may  be  expressed  by 


v=(J^-^T(g-Jttu-J,#), 
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(A3A) 


which,  of  course,  is  contingent  upon  the  invertibility  of  the  square  matrix,  J^J„.  Fur- 
ther, before  estimating  \vj\,  it  is  convenient  to  utilize  the  summing  convention  of  re- 
peated indices,  as  well  as  to  define 

D  =  (J*!.)-1^.  (A3.5) 

Applying  absolute  values  signs  and  the  triangle  inequality,  a  strict  upper  bound  for  \vj\ 
may  be  found  from 

M  <  \Djil\gi\  +  iDjiJ^WSfo]  +  IDjiJ^ikWukl.  (43.6) 

This  representation  assumes  that  all  contributions  to  the  uncertainty  involve  systematic 
errors,  not  random  errors.  An  expression  of  uncertainty  may  be  written  which  assumes 
that  the  contributions  are  treated  as  a  combination  of  random  and  systematic  errors. 
The  total  uncertainty  for  the  ,7th  varied  parameter,  bj,  is  estimated  by 

y/Dji  ('J  +  S1J2J  +  \DjiJ<t>,i\\t>K*\  +  \DjiJu,ik\M,  U3-7) 

where  a  summation  is  implied  on  the  i  index,  and  where  the  variance  in  the  estimated 
deviation  from  the  fit  is 

2  =  ST8  (A3<8) 
9     2M-N  ■  } 

and  the  standard  deviation  in  the  angle  of  incidence  obtained  from  calibration  experi- 
ments of  the  instrument  is  given  as 

s„  =  0.002°.  (A3.9) 

The  systematic  component  of  uncertainty  in  the  angle  of  incidence  is  given  by  8<f>ayt 
which  is  set  equal  to  ±0.002°  from  Table  I.  Equation  (A3. 7)  was  used  to  calculate  the 
relative  contributions  to  the  uncertainties  in  the  model  parameters  for  the  two-layer 
model  reported  in  Table  IV.  The  collective  uncertainty  for  the  sum  of  the  .;'th  and  A;th 
varied  model  parameters,  bj  -f  fet,  is  of  the  form 
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+  (D2ji  +  Dli  +  2DjiDik) 

+  (\DjiJu,ii\  +  \DkiJU)il\)  |u/|  .  (A3.10) 

The  uncertainty  in  the  total  thickness,  if  +  t,-,  as  given  in  Table  IV,  is  calculated  using 
eq.  (A3.10). 
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Appendix  4 


Goodness  of  Fit 

A  statistical  procedure  for  testing  the  goodness  of  fit  of  the  measurement  data  to  a  model 
is  explained  in  reference  [A4.1].  The  procedure  assumes  that  there  are  m  experimental 
determinations  of  A  and     for  estimating  the  n  parameters  of  the  model  and  that  some 
of  these  determinations  represent  multiple  measurements  taken  on  the  same  specimens 
over  time.  The  repeated  measurements  represent  the  true  precision  of  the  process  and 
are  equivalent  to  the  long-term  instrumental  uncertainties  as  explained  in  section  V. 

A  sum  of  squares  (SS)  representing  true  precision,  sometimes  called  a  pure  error  term, 
is  pooled  over  specimens  with  repetitions.  Specifically,  given  mj  repetitions  on  the  ith 
specimen,  the  standard  deviations  associated  with  the  A's  and  -0's,  denoted  by  Si(A) 
and  with  degrees  of  freedom  (mj  —  1)  and  (mj  —  1),  respectively,  are  pooled  with 

standard  deviations  from  the  other  specimens  as  follows: 

SSpure  error  =  (m,  -  1)  (sl(A)  +  slW)  +  (ma  -  1)  (^(A)  +  s\^))  +  . . .  (A4.1) 

Given  the  proper  model,  the  least-squares  fit  to  the  ellipsometric  measurement  data 
will  produce  a  sum  of  squares  for  the  fit  which  is  comparable  to  the  pure  error  sum  of 
squares.  Disagreement  between  the  two  is  taken  as  evidence  that  the  model  does  not  ad- 
equately describe  the  data.  Formally,  this  procedure  partitions  the  total  error  sum  of 
squares  into 

SStotal  error  —  SSpure  error  +  SSlack  of  fit,  (A4.2) 
where  the  SStotai  error  is  defined  in  the  usual  way  as  the  sum  of  squares  of  deviations  or 


SStotal  error  =  £  (AJ  -  A?)2  +  £>J  -        ,  (A4.3) 

*=i 


where  the  measured  values  are  denoted  by  primes  and  the  predicted  values  are  denoted 
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by  double  primes. 

Then,  an  F  statistic  is  constructed  to  test  the  significance  of  the  lack-of-fit  term  as 


where 


jp          \SSt0tai  error       SSpure  error]  /  [df total  error       dfpUre  error]  (.4.4  4) 

[SSpure  error]  /  [dfpure  error] 


dfpure  error  =  2  (mi  -  l)  +  2  (m2  -  l)  + 


and 


DFtotal  error  —  (2m  —  7l). 

The  F  statistic  is  compared  with  Fa(u\,U2)  which  represents  the  upper  a  percent  point 
of  the  F  distribution  with  v\  degrees  of  freedom  in  the  numerator  and  1/2  degrees  of  free- 
dom in  the  denominator.  If 


F<Fa{uuu2),  {A4.5) 

there  is  no  significant  lack  of  fit  at  the  a  significance  level. 

Reference 

A4.1  Draper,  N.,  and  Smith,  H.,  Applied  Regression  Analysis,  second  edition  (John  Wiley 
&  Sons,  New  York,  1981),  pp.  31-38. 
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Appendix  5 

Computer  Program  for  Determining  the  Thickness  of  the  SRM  Silicon  Dioxide  Layer 

A  computer  program  (PROG260),  which  is  a  combination  of  McCrackin's  FORTRAN 
program  [A5.1]  and  a  least-squares  method  for  determining  the  SRM  silicon  dioxide 
layer  thickness,  is  given  on  an  accompanying  5.25-in.  floppy  disk.  McCrackin's  method 
is  used  to  obtain  a  layer  thickness  and  the  least-squares  part  of  the  program  is  then  used 
to  minimize  the  sum  of  the  squares  of  the  deviations  of  A  and  ifr  about  this  thickness. 
The  disk  has  been  formatted  on  an  IBM  PC/XT  and  contains  the  source  code,  object 
file,  and  executable  file  written  using  the  Ryan-McFarland  FORTRAN  Version  2  for  IBM 
PC/XT,  AT,  and  compatible  computers  with  a  math  coprocessor.*  This  program  en- 
ables the  user  to  calculate  the  silicon  dioxide  layer  thickness  using  the  values  of  the  ellip- 
sometric  parameters  A,  if),  and  the  angle  of  incidence,  <f>,  obtained  from  a  single  ellipso- 
metric  measurement. 

The  program  contains  fixed  values  for  the  index  of  the  medium  (air)  and  wavelength, 
A  =  632.8  nm.  It  also  contains  default  values  for  the  angle  of  incidence,  AI  (degrees),  re- 
fractive index  of  the  silicon  substrate,  NS,  refractive  index  of  the  silicon  dioxide  layer, 
NF,  and  an  interlayer  thickness  (nanometers)  and  refractive  index,  INTER.  The  de- 
fault values  are  as  follows:  AI  =  70.0,  NS  =  3.875,  0.018,  NF  =  1.461,  and  INTER  = 
1.0,2.8,0.0.  These  default  values  can  be  modified  according  to  the  user's  requirements. 

For  the  two-layer  model,  the  determination  of  the  thickness  of  the  silicon  dioxide  layer  is 
based  upon  using  the  best  estimates  of  the  refractive  index  of  the  silicon  dioxide  layer, 
refractive  index  and  thickness  of  the  interlayer,  and  refractive  index  of  the  substrate. 
The  procedure  is  as  follows.  The  two-layer  model  data  from  the  certificate  are  entered 
with  the  data  entry  commands,  AI,  NF,  INTER,  and  NS.  Entering  the  principal  angle, 
^,  using  the  AI  command  and  entering  the  values  for  A  and     using  the  CD  command 
will  produce  the  thickness  of  the  silicon  dioxide  layer.  By  adding  the  value  reported  as 
the  interlayer  thickness,  1.0  nm,  as  it  is  reported  on  the  certificate,  to  the  calculated 
thickness,  if  +  t{  will  be  obtained.  The  calculation  for  the  one-layer  model  is  performed 
in  the  same  manner,  but  with  the  thickness  and  refractive  index  for  the  interlayer  set  to 
zero  by  entering  the  command  INTER  and  the  data  "0.0."  The  one-layer  refractive  in- 


dex  must  be  entered  from  the  certificate  using  NF  and  the  appropriate  value  (this  infor- 
mation can  be  found  on  page  3  of  the  certificate).  The  CD  command  will  then  produce 
the  layer  thickness  for  the  one-layer  model. 

Reference 

A5.1  McCrackin,  F.L.,  A  FORTRAN  program  for  analysis  of  ellipsometer  measurements, 
NBS  Technical  Note  479  (1969). 
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